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ABSTRACT 

We present new models for the formation of disk galaxies that improve upon previ- 
ous models by following the detailed accretion and cooling of the baryonic mass, and 
by using realistic distributions of specific angular momentum. Under the assumption 
of detailed angular momentum conservation the disks that form have density distri- 
butions that are more centrally concentrated than an exponential. We examine the 
influence of star formation, bulge formation, and feedback on the outcome of the sur- 
face brightness distributions of the stars. Low angular momentum haloes yield disk 
galaxies with a significant bulge component and with a stellar disk that is close to 
exponential, in good agreement with observations. High angular momentum haloes, 
on the other hand, produce stellar disks that are much more concentrated than an ex- 
ponential, in clear conflict with observations. At large radii, the models reveal distinct 
truncation radii in both the stars and the cold gas. The stellar truncation radii owe to 
our implementation of star formation threshold densities, and are in excellent agree- 
ment with observations. The truncation radii in the density distribution of the cold 
gas reflect the maximum specific angular momentum of the gas that has cooled. We 
find that these truncation radii occur at HI surface densities of roughly 1 M© pc~^, 
in conflict with observations. We examine various modifications to our models, in- 
cluding feedback, viscosity, and dark matter haloes with constant density cores, but 
show that the models consistently fail to produce bulge-less disks with exponential 
surface brightness profiles. This signals a new problem for the standard model of disk 
formation: if the baryonic component of the protogalaxies out of which disk galaxies 
form have the same angular momentum distribution as the dark matter, disks are too 
compact. 

Key words: galaxies: formation — galaxies: fundamental parameters — galaxies: 
spiral — galaxies: kinematics and dynamics — galaxies: structure — dark matter. 



1 INTRODUCTION 

In the current paradigm for galaxy formation, set forth by 
White & Rees (1978) and Fall & Efstathiou (1980), galaxies 
are considered to form through the cooling of baryons inside 
dark matter haloes that grow by means of gravitational in- 
stability and acquire angular momentum from cosmological 
torques. If the cooling baryons conserve their specific angu- 
lar momentum, disk galaxies will form in the centres of the 
haloes, with scale lengths that are in good agreement with 
observations (e.g., Dalcanton, Spergel & Summers 1997; Mo, 
Mao & White 1998; de Jong & Lacey 2000). 

However, detailed hydrodynamical simulations aimed 
at investigating this process of galaxy formation have indi- 
cated an important problem for the cold dark matter (CDM) 
model. In CDM cosmologies haloes form hierarchically by 
the merging of many lower mass haloes. Because the cooling 
in dense, low mass haloes is very efficient, the baryons in 



these systems have already cooled by the time they merge 
with the more massive protogalaxy. They reach the centre of 
the potential well by means of dynamical friction, through 
which they loose a significant fraction of their specific angu- 
lar momentum to the dark matter. Consequently, the disks 
that form in these simulations are an order of magnitude 
too smaU (Navarro & Benz 1991; White & Navarro 1993; 
Navarro & Steinmetz 1999). This has become known as the 
angular momentum catastrophe of disk galaxy formation. 



Several solutions to this problem have been suggested. 
Weil, Eke & Efstathiou (1998), Domi'nguez-Tenreiro, Tissera 
& Saiz (1998), Sommer-Larsen, Gelato & Vedel (1999) and 
Eke, Efstathiou & Wright (2000) argued that stellar feed- 
back and/or ionizing background radiation can prevent the 
cooling of gas in these small mass haloes, therewith consider- 
ably reducing the angular momentum loss (but see Navarro 
& Steinmetz 1997) . An alternative suggestion has been to al- 
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ter the power spectrum of initial density fluctuations, either 
by invoking an alternative form of dark matter (Sommer- 
Larsen & Dolgov 2001), or by resorting to a specific model 
for inflation (Kamionkowski & Liddle 2000). At the present, 
it is still unclear which of these solutions, if any, is most 
successful. What has become clear, however, is that if cos- 
mological torques are indeed the main source of angular mo- 
mentum on galactic scales, there is very little room for any 
angular momentum loss if we want to explain the sizes of 
present day disk galaxies. 

Motivated by this understanding, the past couple of 
years a large number of studies have presented analyti- 
cal models for the formation of disk galaxies that rely on 
the assumption of detailed angular momentum conserva- 
tion. Kauffmann (1996) investigated the properties of disk 
galaxies within this framework, and linked it to the evo- 
lution of damped Lya absorption systems. Dalcanton et 
al. (1997) and Mo et al. (1998) investigated the structural 
properties of disks, with emphasis on the variance induced 
by the distribution in halo angular momentum. Subsequent 
studies included recipes for bulge formation, gas viscosity, 
star formation and/or feedback and investigated more de- 
tailed properties of these model disk galaxies, such as the 
TuUy-Fisher relation, the gas mass fractions, and the ori- 
gin of the Hubble sequence (van den Bosch 1998, 2000; 
Jimenez et al. 1998; Natarajan 1999; Heavens & Jimenez 
1999; van den Bosch & Dalcanton 2000; Firmani & Avila- 
Reese 2000; Avila-Reese & Firmani 2000; Efstathiou 2000; 
Zhang & Wyse 2000; Buchalter, Jimenez & Kamionkowski 
2001; Ferguson & Clarke 2001). 

One important shortcoming of most of these models, 
however, is that they make the a priori assumption that the 
cooled gas arranges itself in an exponential disk. However, 
one of the open issues in the formation of disk galaxies, is 
to actually understand why they reveal a universal surface 
brightness distribution. If disk galaxies indeed form as envi- 
sioned in our standard picture, their resulting density distri- 
bution is directly related to the specific angular momentum 
distribution of the protogalaxy. Motivated by the work of 
Mestel (1963) and Crampin & Hoyle (1964), Dalcanton et 
al. (1997) made the assumption that the protogalaxy has 
the angular momentum distribution of a uniform sphere in 
solid-body rotation, and showed that the resulting disks are 
actually more centrally concentrated than an exponential. 
Firmani & Avila-Reese (2000) used more realistic distribu- 
tions of mass and angular momentum for the dark matter 
haloes, but again found similarly concentrated disks. Bullock 
et al. (2000) determined the distribution of specific angular 
momentum in CDM haloes from high resolution TV-body 
simulations, and again concluded that these distributions 
will form overly concentrated disks. Therefore, if our picture 
for disk formation is correct, the question arises whether the 
subsequent processes of star formation, bulge formation and 
feedback can produce exponential stellar disks in agreement 
with observations. 

Here we present new, but similar, models for the forma- 
tion of disk galaxies. These models, presented in Section H 
will be used in forthcoming papers to investigate a wide 
range of disk properties such as colors and metallicities, gas 
mass fractions, and the TuUy-Fisher relation. In this paper, 
however, we focus solely on the density distributions of the 
resulting disks. In Section H we investigate the effects of star 



formation and bulge formation on the outcome of the surface 
brightness profiles of the disks. In Section MJ we compare the 
truncation radii predicted by the models with observations. 
A problem with reproducing low surface brightness galaxies 
is discussed in Section H, and we summarize our results in 
Section ra. 



2 THE MODELS 

2.1 The basic framework 

The main assumptions that characterize the framework of 
our models are the following: (i) dark matter haloes around 
disk galaxies grow by the smooth accretion of mass, (ii) the 
angular momentum of protogalaxies originates from cosmo- 
logical torques, (iii) in the absence of cooling the baryons 
have the same distribution of mass and angular momentum 
as the dark matter, and (iv) the cooling baryons conserve 
their specific angular momentum. 

The two main ingredients that determine the formation 
and evolution of a disk galaxy, therefore, are (the evolu- 
tion of) the mass and angular momentum of the virialized 
object; Mvir(r, «) and Jvir(r, z)n. We characterize the an- 
gular momentum of the protogalaxies by the dimensionless 
spin parameter A — Jvir|-Evir|'^''^G~^MJT^' . Here -Evir is the 
halo's energy, and G is the gravitational constant. We follow 
Firmani & Avila-Reese (2000) and make the additional as- 
sumptions that (i) the spin parameter A of a given galaxy is 
constant with time, (ii) each mass shell that virializes is in 
solid body rotation, and (iii) the rotation axes of all shells 
are aligned. 

None of these assumptions are necessarily accurate. For 
instance, the assumption of smooth mass accretion seems 
inconsistent with the hierarchical merger picture of struc- 
ture formation in a CDM universe. On the other hand, the 
fragility of disks (e.g., Toth & Ostriker 1992) suggests that 
mergers can not have played a dominant role in establishing 
the main properties of disk galaxies. Furthermore, numerical 
simulations suggest that too much merging results in disks 
that are too small, as baryons tend to loose their angular 
momentum to the dark matter in the process. The assump- 
tions regarding the angular momentum are also question- 
able. Nevertheless, as we show in Section 5.5, these assump- 



tions result in haloes with an angular momentum profile that 
is in good agreement with high resolution A''-body simula- 
tions. Finally, we like to stress that the main goal of this 
study is to explore how the final disk characteristics depend 
on the various model ingredients. We therefore adhere to 
simple, parameterized descriptions of the mass accretion his- 
tory (hereafter MAH) and angular momentum distribution. 
Although perhaps not completely realistic, they provide use- 
ful insights. 

The main outline of the models is as follows. We set 
up a radial grid between r = G and the present day virial 
radius of the model galaxy and we follow the formation and 
evolution of the disk galaxy using a few hundred time steps. 
We consider six mass components: dark matter, hot gas, disk 
mass (both in stars and in cold gas), bulge mass, and mass 



* Throughout this paper, r and z refer to spherical radius and 
redshift, respectively. 
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ejected by outflows from the disk. The dark matter, hot gas, 
and bulge mass are assumed to be distributed in spherical 
shells, whereas the disk stars and cold gas are assumed to 
be in infinitesimally thin annuli. Each time step we compute 
the changes in these various mass components in each radial 
bin. Below we describe the detailed prescriptions used. 

2.2 The evolution of the dark matter component 

The backbone of the models is the formation and evolution 
of the dark matter haloes, which is determined by the pa- 
rameters of the background cosmological model and by the 
power spectrum P{k) of the initial density fluctuations. The 
parameters of importance for the models presented here are 
the present day matter density flo, the present day density 
in the form of a cosmological constant Qa, the Hubble con- 
stant h = 1/0/(100 km s~^ Mpc~^), the baryon density fibar, 
and the normalization erg of P{k). 

As discussed above, we make the assumption that proto- 
galaxies accrete mass smoothly. Rather than attempting to 
link the actual accretion rate to the cosmological framework, 
using for instance the extended Press-Schechter formalism 
to construct merger histories, we adopt a simple parame- 
terization. This has the advantage that we can describe the 
mass accretion history (hereafter MAH) by one or two free 
parameters. 

For a given virialized mass at 2 = 0, Mvir(O), we write 
the MAH as 

-1 l/a 



M^i,{z) = Afvir(O) 



1 



log(l 



log(l + Zf) 



(1) 



Here a and 2/ are free parameters. The parameter a de- 
scribes whether most of the flnal virial mass is accreted 
early or late. In Figure YU we plot MAHs for zj — 10 and 
three different values of a: 0.3, 0.6 and 0.9. The correspond- 
ing redshifts at which half the present day mass is assem- 
bled are 0.57, 1.26, and 2.04, respectively. A comparison 
with MAHs computed using the extended Press-Schechter 
approximation based on the conditional probabilities for a 
Gaussian random field shows that the parameterization (M 
with Zf = 10 and 0.3 ^ a ^ 0.9 comprises ~ 95 percent 
of the typical MAHs for a typical galaxy sized halo, with 
a = 0.6 close to the average (cf. Lacey & Cole 1993; Fir- 
mani & Avila- Reese 2000; Buchalter et al. 2001). Clearly, in 
reality the values of 2/ and a will depend on the mass of the 
object. This is not taken into account here. The aim of this 
paper is not to use the most accurate MAHs, but merely to 
use a simple parameterization to investigate how the result- 
ing disk surface densities depend on the MAH. Therefore, in 
what follows, we use values of a in the range 0.3-0.9 (with 
Zf fixed at 10) to gauge the dependence of our models on 
the MAHs. When required, we shall explore a wider range 
of Zf and a values. 

The virial radius of dark matter haloes as function of 
redshift is given by 

A/vir(2) V''Yt2(2)^'^' 
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Figure 1. The mass accretion histories (MAHs) defined by equa- 
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for z 



f 



10 and three different values of a. These curves 



describe the growth of halo mass, normalized to the present day 
mass, as function of redshift. A comparison with Figure 1 in Fir- 
mani & Avila-Reese (2000) shows that these MAHs are in reason- 
able agreement with the 5, 50, and 95 percentile MAHs computed 
using the extended Press-Schechter approximation. 



inside r^h expressed in terms of the critical density for clo- 
sure. We use the fitting formula of Bryan & Norman (1998): 



Avir(2) = 187r^ + 82a; - 39a; 



(3) 



with X = fl{z) — 1. This fitting formula is valid for a fiat 
Universe (f7o + ^a = 1) and accurate to one percent in the 
range 0.1 < n{z) < 1.0. 

Rather than attempting to compute the density dis- 
tribution of virialized haloes from extended secondary infall 
models, as for instance in Avila-Reese, Firmani & Hernandez 
(1998), we adopt a simple parameterized density distribu- 
tion. We assume that the dark matter virializes such that at 
each redshift the halo's density distribution is given by 



Pv.(r)=p.(^)"^(l + ^) 



7-3 



(4) 



with Avir the virial density, defined as the average density 



with ps and r^ dependent on 2 and Mviriz). This profile 
is motivated by both numerical simulations and by the ob- 
served rotation curves of disk galaxies: numerical simula- 
tions of structure formation in a CDM universe suggest val- 
ues for 7 in the range 1.0 — 1.5 (e.g., Navarro, Frenk & White 
1996, 1997; Fukushige & Makino 1997; Moore et al. 1998), 
whereas rotation curves of disk galaxies imply < 7 jj 1.5 
(van den Bosch et al. 2000; van den Bosch & Swaters 2001). 
We mainly focus on models with 7 = 1 for which equa- 
tion (W reduces to the NFW proflle (Navarro, Frenk & 
White 1997). However, because of the relatively large un- 
certainties we shall also consider other values of 7. 

The total mass, energy, and angular momentum of a 
halo with density distribution (M) are given by: 

3 J 



M 



47rpsr^/(c) 



(5) 
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Svir = -iMvirV;i/(c)"'/l(c), (6) 

and 

Jvir = V2 ArvirMvirKir/(c)/l(c)"'''^ (7) 

Here Vvir is the circular velocity of the halo at the virial 
radius, c = rvii/vs is the halo concentration parameter, and 
f{x) and h{x) are given by 



/(x)= / dyy^-^il + yy-^ 
Jo 

and 



h{x) 



dyf{y)y'-^{l + y) 



N7-3 



(8) 



(9) 



Equations ([lJ)-(0) completely specify the dark matter den- 
sity distribution at each redshift once we know the value of 
the concentration parameter c as function of z. We use the 
model of Bullock et al. (2001), which gives 



c{M^ 



= Co 



/ 1 + Zcoll \ 
\ l + Z J 



(10) 



with ZcoU the redshift at which a mass / • Mvir collapses. 
Here Co and / are constants that depend on the cosmology. 

Note that in our approach all haloes at the same red- 
shift and with the same mass have the same concentration 
parameter c. However, in reality there is a large spread in 
values of c (e.g., Jing 2000; Bullock et al. 2001). Furthermore 
one expects that c is correlated with the actual MAHs, such 
that MAHs with a higher value of a (i.e., more early accre- 
tion) are more concentrated (e.g., Avila- Reese, Firmani & 
Hernandez 1998; Firmani & Avila- Reese 2000). We ignore 
this scatter in c and its associated correlation with a. De- 
spite these obvious shortcomings, the simplicity of the above 
parameterizations allows us to investigate the main depen- 
dencies of our models, which is our main concern in this 
paper. 

Numerous studies, both analytical and numerical, have 
shown that the distribution of the spin parameter A is well 
described by a log-normal distribution 



p{\)d\ ■ 



1 



crAv27r 



exp 



ln2(A/A) 



2al 



dA 
A ' 



(11) 



(e.g., Barnes & Efstathiou 1987; Ryden 1988; Cole & Lacey 
1996; Warren et al. 1992). In this paper we gauge the A- 
dependence of our models by constructing models with A = 
0.028, 0.06, and 0.129. These values correspond to the 10, 50, 
and 90 percentile points of the distribution of equation ( |ll[ ) 
with A = 0.06 and ax ~ 0.6. 



2.3 Disk formation 

In order to compute the formation and evolution of the disks 
we proceed as follows. Each time step At a shell with mass 
AM = Mvir(i) — Mvir(i — At) virializes. A fraction /bar = 
fibar/f^o of this mass is in baryons, and is heated to the 
halo's virial temperature 



_ 1 /imp 2 



(12) 



where fj,mp is the mass per particle, and k is Boltzmann's 
constant. This gas is added to the hot gas component be- 
tween radii rmin and rmax, whereby we assume that the gas 



follows the same density distribution as the dark matter. 
Assuming no shell crossing we add AAf to the outer parts 
of the halo, and we thus set rmax equal to the virial radius 
at t, and rmin to the radius inside of which the total mass 
(baryons plus dark matter) is equal to the virial mass at 
time t — At. The baryons dissipate energy radiatively, but 
are assumed to conserve their specific angular momentum. 
The time scale on which they will reach centrifugal equilib- 
rium in the disk is given by tc = max[iff , tcooi]- Here tg is 
the free-fall time defined as 



tff = 



32Gp' 



(13) 



(14) 



with p the average halo density, and 

, O Klv'iY pe 

f-cool — 77 /^^^P 7 7"^ T T 

2 pgasAjv^Zhotj pe — i 

is the cooling time. Here fie is the number of particles per 
electron, and Ajv(.^hot) is the normalized cooling function 
for a gas with metallicity Zi^^t ■ For Ajv we use the coUisional 
ionization equilibrium cooling functions of Sutherland & Do- 
pita (1993), assuming a Helium mass abundanc e of 0.25. 
Within our framework outlined in Section 2.1, the ra- 



dius at which the baryons reach centrifugal equilibrium is 
computed as follows. For a shell with rmin < r < rmax in 
solid body rotation with circular frequency uo and with den- 
sity distribution p(r) the total angular momentum is given 
by 



Jslicll = — TT OJo 



dr r p(r) 



(15) 



One can also write the angular momentum of the shell as 
AJ = Jvir (t) — Jvir (i — Ai), 1.6., the difference in the total an- 
gular momentum of the halo before and after the addition of 
the new shell of matter. Here Jvir(t) is given by equation {m. 
Thus, for the shell's circular frequency one finds 



Wo 



_3_ 

87r 



AJ 



dr r p(r) 



(16) 



At time t' = t + tc we add these baryons to the disk. Us- 
ing detailed conservation of specific angular momentum, we 
compute the baryonic mass added to the disk annulus with 
rk < r < rfe+i as the baryonic shell mass with cylindrical 
radii Rk < R < Rk+i- Here 



Ri 



riVc(ri,t') 

LOO 



(17) 



with Vc{ri,t') the total circular velocity of the galaxy at 
radius ri at time t' . Throughout we assume that the disk is 
infinitesimally thin, and each time step we use the adiabatic 
invariant formalism of Blumenthal et al. (1986) and Flores 
et al. (1993) to compute the gravitational contraction of the 
dark matter induced by the baryons settling in the disk. 

2.4 Star formation 

Star formation in (disk) galaxies is a complicated process 
which is only poorly understood. In particular, it is unclear 
whether, in isolated disk galaxies, star formation is mainly 
triggered by density waves (Wyse 1986; Wyse & Silk 1989), 
by local feedback from supernovae (Gerola & Seiden 1978), 
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or by cloud-cloud interactions (Tan 2000). Furthermore, al- 
though it is generally accepted that the star formation is 
self-regulated, the detailed interplay between gravitational 
instabilities, supernovae feedback, turbulent viscosity, mag- 
netic fields, and star formation is still heavily debated (see 
e.g., Firmani & Tutukov 1994; Silk 1997; Tan 2000; and ref- 
erences therein). Therefore, we use a star formation recipe 
that is motivated by empirical findings, so that despite our 
ignorance, we are implicitly taking account of all these de- 
tailed physical processes. 

Observations have shown that the star formation rates 
in disk galaxies are well fit by a simple Schmidt (1959) law: 



E = risN AM, EsN 



(23) 



Tp = £SF Sg 



(18) 



This simple empirical law holds over many orders of mag- 
nitude in gas surface density, and even applies to circum- 
nuclear starburst regions. However, when applied to local 
gas densities, the Schmidt law breaks down at large disk 
radii, where the star formation is found to be abruptly sup- 
pressed. In a seminal paper, Kennicutt (1989) showed that 
these radii correspond to the radii where the gas surface den- 
sity falls below the critical surface density given by Toomre's 
(1964) stability criterion: 

Tgas k{R) 



Scrit(-R) 



(19) 



3.36 GQ 

Here Q is a dimensionless constant near unity, (Tgas is the 
velocity dispersion of the gas, and k is the epicycle frequency 
given by 

1/2 



k(R) = V2 



VcjR) 
R 



1 + 



R dVc 
Vc{R) dR 



with Vc{r) the circular velocity of the entire system. 
Solving dEgas/di — — f/), yields 



Egas(i) — 



m 



At) 



(20) 



(21) 



with m = 1/(1 — n) (cf. Heavens & Jimenez 1999). In each 
annulus in the disk we then compute the mass in stars 
formed between t — At and t as 

AM. = A [Egas(t - At) - Sgas(f)] , (22) 

with A the area of the annulus and Egas(t) = 
max[Ecrit(t), Egas(t)]. This way, star formation is not al- 
lowed to deplete gas to surface densities below Ecrit. 



2.5 Feedback by Supernovae 

When stars evolve they put energy into the interstellar 
medium (ISM) which impacts on the further evolution of the 
galaxy. By resorting to an empirical description of the star 
formation, we are implicitly taking account of the eff'ects 
that these feedback processes have on the star formation 
rate. What is not taken into account, however, is a possible 
feedback-driven outfiow of gas from the disk. Here we use a 
simple parametric model, similar to the ones used in vari- 
ous semi-analytical models for galaxy formation. We assume 
that the amount of gas blown out of the disk is proportional 
to the total energy input by supernovae (SNe) and inversely 
proportional to the escape velocity squared. At each time 
step we compute the total energy injected into the ISM at 
each disk annulus as 



Here AM* is the mass in stars formed (equation |22[|), 
EsN = 10^^ ergs is the energy produced by one SN, and 
»7SN is the number of SNe per solar mass of stars formed. 
We assume that this energy drives a galactic wind, whereby 
a mass AMejoct is blown out of the halo. Requiring energy 
balance one obtains 

2 £fb »7SN EsN 



AMei 



■AM. 



(24) 



(cf. Kauffmann, White & Guiderdoni 1993; Natarajan 1999). 
Here Vise is the local escape velocity, and £fb is a free pa- 
rameter that describes what fraction of the energy released 
by SNe is converted into kinetic energy to drive the outfiow. 
For simplicity, we assume that the ejected mass is forever 
lost from the system: the ejected mass is not considered for 
later infall, and the corresponding metals are not used to 
enrich the infalling gas. 

2.6 Bulge formation 

Self-gravitating disks tend to be unstable against global in- 
stabilities such as bar formation. Here we follow the ap- 
proach of van den Bosch (1998, 2000) and Avila-Reese & 
Firmani (2000) and assume that an unstable disk trans- 
forms part of its disk material into a bulge component in a 
self-regulating fashion such that the final disk is marginally 
stable. 

Motivated by the work of Christodoulou, Shlosman & 
Tohline (1995) we consider a disk to be unstable if 

/Vdisk(j-)\ ^ fr,r\ 

Qmax = max — — --1- < acrit. (25) 

0<r<r^i, \Vciyc{r) J 

Here Ofcrit is a free parameter, which regulates the disk-to- 
bulge ratios of the final model galaxies, and Vdisk(T') and 
Vcirc('") are the circular velocities of the disk (cold gas plus 
stars) and the composite disk-bulge-halo system, respec- 
tively. 

The formation of bulges out of unstable disk material 
is a complicated process. It is likely to involve bars, which 
are efficient in transporting gas inwards, and which subse- 
quently dissolve to form the bulge. However, the details of 
the mass fiow, and the shape parameters of the resulting 
bulge component are poorly understood, and we are forced 
to make some ad hoc assumptions. If the disk is unstable 
according to (Eq), we transform cold gas mass from the in- 
side out to the bulge, whereby we assume that the infiow is 
such that it does not create a positive gradient in the surface 
density of the cold gas. This particular choice for extracting 
bulge material from the disk is based on extensive tests, and 
is optimized to yield stellar disks that are close to exponen- 
tial. We s hort ly address the infiuence of these assumptions 
in Section 3.4. 



At each time step we use an iterative procedure to com- 
pute the gas mass transformed to the bulge, AMbuige, such 
that the resulting disk has amax = ctcrit . This mass AMbuigc 
is assumed to form stars instantaneously with 100 percent 
efficiency, and the SN energy released by this burst of star 
formation is added to the energy released by the quiescent 
star formation in the disk at the disk radii from which the 
bulge material originates. This latter assumption has no im- 
portant consequences for our results. If, for instance, we were 
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to deposit all SN energy related with the bulge formation at 
r = 0, this does not alter any of our conclusions. 

We model the bulge as a sphere with a Hernquist den- 
sity profile: 



Pbi^) = ^ 



Tb 



27r r {r + rt)^ ' 



(26) 



where rt is a scale length (Hernquist 1990). Given our poor 
understanding of the detailed processes involved we use em- 
pirical relations to compute r;,. Andredakis, Peletier & Bal- 
cells (1995) have shown that the effective radius Ve (defined 
as the radius encircling half of the projected light) is directly 
related to the total B-band luminosity of the bulge by the 
empirical relation 



Mb = -19.75- 2.8 log(re). 



(27) 



We use this relation to compute the scale length of the 
bulge, taking account of the fact that for a Hernquist sphere 
re = 1.8153 r;, (Hernquist 1990). In practice, our results do 
not depend on this simple scaling assumption. The main pa- 
rameter for the bulge is its total mass; changes in its actual 
density distribution are only a second-order effect, and do 
not influence our results. 



2.7 Stellar population modeling & chemical 
evolution 

In order to convert the stellar masses into luminosities we 
use the latest version of the Bruzual & Chariot (1993) stel- 
lar population synthesis models. These models provide the 
luminosities li{t, Z) of a single burst stellar population with 
a total mass of 1 Mq as function of age t and metallicity 
Z in various optical passbands i. In each disk annulus we 
keep track of the amount of stars formed at each time step, 
AM«(tfe), as well as the metallicity of the cold gas Z{tk)- 
This allows us to compute the total luminosity in that an- 
nulus, at any time step tj, as 



Li [tj 



3 

E 



h[tj-tk,Z{tk)]/^M,{tk) 



(28) 



In order to model the chemical enrichment of the ISM 
we follow the standard instantaneous recycling approxima- 
tion (IRA). We assume that a fraction TZ of the mass in stars 
formed is instantaneously returned to the cold gas phase 
with a yield y (which is defined as the fraction of mass con- 
verted into stars that is returned to the ISM in the form of 
newly produced metals). In each disk annulus, and at each 
time step, mass conservation implies 



AMcoid = AAf.ooi - (1 - 7^)AA/, - AA/ojoct 
and for the mass in metals one thus obtains 



AAfn 



= Z-hot^Mc, 



iAAf, 



eject 



^cold(l-7^)AA/. +j/AAf, 



(29) 



(30) 



We use these two equations to track the evolution of the 
metallicity of the cold gas in the disk, Z^oid , as function of 
both time and radius. Note that the mass that is ejected is 
forever lost from the galaxy. The associated metals are not 
used to enrich the hot gas in the halo, and Zhot is therefore 
constant with time. 



2.8 Detailed description of model parameters 

We distinguish two different sets of parameters. The first set, 
which we call the galaxy parameters, consists of Afvir(O), A, 
and a, and specify a particular model galaxy. The second set, 
the model param,eters, consists of parameters that specify 
the particular formation model. The model parameters are 
as follows: 

• Cosmology: f7o, f^A, f^bar, h, as. In this paper we re- 
strict ourselves to the currently popular ACDM cosmology 
with 9,0 = 0.3, f^A = 0.7, h = 0.7, and as = 1.0. These 
parameters are currently favored by a large body of obser- 
vations. For the baryon density we adopt fibar = 0.019 h~^ 
as suggested by the observations of primordial deuterium 
abundances by Tytler et al. (1999). 

• Halo structure: Co, /, and 7. For our ACDM cosmology, 
Bullock et al. (2001) found that dark matter haloes in N- 
body simulations have 7 = 1.0, Co = 4.0 and / = 0.01. 
Unless stated otherwise, we shall adopt these values. 

• Star formation: esF, n, Q, and cTgas. For a large sample 
of disk galaxies and circum-nuclear starburst regions Ken- 
nicutt (1998) found esF = 2.5 x lO""* Moyr'^kpc"^ and 
n = 1.4, while it was shown in Kennicutt (1989) that for 
fgas = 6 kms~^ and Q = 1.5 the star formation truncation 
radii are in good agreement with observations. Since our star 
formation recipe is based on empirical relations, we keep our 
star formation parameters fixed at these observationally de- 
termined values. 

• Feedback: £fb. Since there are no empirical constraints 
on this feedback efficiency parameter we consider efb a free 
parameter. 

• Bulge formation: acrit. This parameter sets the amount 
of self-gravity of the disk above which we transfer cold 
gas from the disk to the bulge. Numerical simulations by 
Christodoulou et al. (1995) found cicrit ~ 0.7 for a gaseous 
disk, which is what we adopt throughout. 

• Stellar populations and chemical enrichment: t;sn, Tl, 
y, Zhot, and the choice for an initial mass function (IMF). 
Throughout we adopt the Scalo (1986) IMF, for which 
7?SN = 4 X 10"^ Mq"^ and 7^ = 0.25. We keep the stel- 
lar yield fixed at y = 0.02, which is a typical value used in 
chemical evolution models. Unless stated otherwise we adopt 
Zhot = 0.3^0, typical of the hot gas in clusters (Mushotzsky 
& Loewenstein 1997). 



3 THE DENSITY DISTRIBUTION OF DISK 
GALAXIES 

3.1 Cooling properties 

We start our investigation into the structural properties of 
our model galaxies by computing the density distributions 
of the cold gas disks in the absence of star formation, feed- 
back, and bulge formation. To that extent we first examine 
the disk mass fractions as function of the various input pa- 
rameters. 

The free-fall time at each redshift z is the same for 
all galaxies as it depends only on the average halo density 
p{z) = Avir(2) Pcrit(2). It iucreases from 0.08 Gyr at 2 = 
10 to 2.2 Gyr at 2 = 0. The cooling time, on the other 
hand, depends on the density and (virial) temperature of 
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Figure 2. The fraction, /coob of baryonic mass inside the virial radius that has cooled and settled in present day disks as function of 
the present day virial mass. Results are shown for three different MAHs, and for two different metallicities of the gas. For low mass 
systems /cool depends only on the free-fall time and is independent of halo mass. For more massive systems /cool is dominated by the 
cooling time scale and /cool decreases rapidly with increasing mass. Accreting more mass at earlier times (i.e., larger a) and increasing 
the metallicity of the gas both result in larger values of /cool ■ 



the gas, and is a strong function of Mvir- In Figure H we plot 
the fraction /cool of baryons inside the virial radius that 
has settled in the disk at z; = as function of the present 
day virial mass. Results are shown for three different values 
of a and for both zero-metallicity gas (left panel) and gas 
with a metallicity one-third Solar (right panel). For gas with 
Zhot = and Mvir(O) ^ lO^/i'^M© one finds that fcooi < iff 
at all redshifts, and between 81 and 93 percent of the gas 
has settled in a disk hy z — 0, depending on the MAH. For 
more massive systems the cooling time exceeds the free-fall 
time during the later stages of evolution, and less mass can 
cool. Consequently, /cool decreases strongly with mass: for 
a system with Mvir(O) — 10^"^h~^ Mq only between 14 and 
44 percent of the baryons inside the virial radius can cool. 
If the gas is pre-enriched to Zhot ~ 0.3^0, the cooling is 
more efficient, and /cool is significantly larger for high mass 
systems. To summarize, /cool is a strong function of mass, 
MAH and metallicity, and one thus expects strong variations 
of the mass fractions of cold gas and stars in galaxies. 



3.2 The density distribution of the gas 

In Figure H we plot the surface density of the cooled gas and 
the total circular velocity at z = as function of radius (no 
star formation or bulge formation is included here). The dif- 
ferent panels show the results of varying one parameter with 
respect to a fiducial model with Mvir(O) = 5 x 10^^h~^ Mq, 
A = 0.06, a = 0.6, and 7 = 1.0 (this model is plotted as 
a solid line in all panels). In the panels on the left we vary 
the spin parameter A. As expected, systems with more an- 
gular momentum (larger A) result in more extended disks. 
For A ^ 0.05 the compact disks that form are strongly self- 
gravitating, resulting in strongly declining rotation curves. 



Once bulge formation is included, these systems will trans- 
form significant fractions of their disk mass in a bulge com- 
ponent, such that the final rota tion curves are in agreement 
with observations (see Section 3.4). The middle panels of 
Figured plot the results for three different MAHs. An earlier 
MAH (larger a) yields disks that are more centrally concen- 
trated, but also more extended. This owes to the differences 
in the cooling histories of the gas. Galaxies with an earlier 
MAH have more gas cooled at high redshifts when the halo 
is denser. Consequently, that gas reaches centrifugal equilib- 
rium at smaller radii. In addition, the final value of /cool is 
higher (cf. Figure 0) , and this excess gas cools to large radii, 
resulting in a more extended disk. Finally, in the panels on 
the right we plot results for three different values of 7. Less 
strongly cusped haloes (lower 7) result in more extended 
disks and lower values of K:irc. This owes to the fact that 
haloes with lower 7 are less dense, such that the gas reaches 
centrifugal equilibrium at larger radii. 

As is immediately apparent from Figure H, none of the 
disk surface density profiles resemble an exponential. In all 
cases the gas is much more centrally concentrated, and is 
well-fitted by a power-law Egas(i?) oc R~" with 1 ^ n ^ 2. 
This confirms previous results by Kauffmann (1996), Dal- 
canton et al. (1997) and Firmani & Avila-Reese (2000), and 
is a direct consequence of the detailed conservation of an- 
gular momentum, and the fact that the dark matter haloes 
have (broken) power-law density distributions (see discus- 
sions in Seiden, Schulman & Elmegreen 1984 and Yoshii 
& Sommer-Larsen 1989). In order to produce stellar disks 
with an exponential surface density distribution out of these 
power-law gas disks one needs to either (i) have specifically 
tuned star formation efficiencies, (ii) transform part of the 
central gas disk into a bulge component, or (iii) somehow 
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Figure 3. The upper panels plot the surface density distribution of cold gas as function of radius for models without star formation and 
bulge formation. The lower panels plot the corresponding circular velocity curves, out to the cut-off radius of the disk. All models have 
the same virial mass of Mvir(O) = 5 X 10^^h~^ Mq. The solid lines correspond to the model with A = 0.06 (median spin parameter), 
a = 0.6 ('average MAH') and 7 = 1.0 (NFW halo profile). The dashed and dotted lines show the influence of varying one parameter as 
indicated in the top panels (see text for a discussion). 



remove the central excess of gas from the disk. Below we ex- 
amine how our particular recipes for star formation, bulge 
formation, and feedback fair in this respect. 

3.3 The influence of star formation 

We now include star formation in o ur m odels using the em- 
pirical parameters listed in Section 2.8 In Figure Wlwe plot 
the present day surface densities of both the gas and stars 
as functions of the normalized radius R/R^i,-. Results are 
shown for six model galaxies that only differ in total mass 
and angular momentum as indicated in the various panels. 

Except for the case with Mvir(O) = 5 x 10^ /i~^ M0 and 
A = 0.129 (lower right corner) all stellar disks are more con- 
centrated than an exponential. Although the density distri- 
butions are reasonably close to exponential over their in- 
termediate radial range, they harbor a pronounced central 
cusp and a distinct cut-off radius. The surface density of the 
cold gas follows Ecrit (indicated by the thin dotted lines), 
except in the outer regions of the disk, where the gas has 
only recently been accreted, and star formation has not had 
sufficient time to deplete the density of the gas to Ecrit . 

A comparison between the upper and lower panels of 
Figure y reveals that the lower mass systems produce disks of 
lower surface brightness and with higher gas mass fractions. 
As pointed out by previous studies, this scaling is in good 



agreement with observations (e.g., Dalcanton et al. 1997; Mo 
et al. 1998; van den Bosch 2000; van den Bosch & Dalcanton 
2000). From Figure it is also apparent that lower mass 
systems have larger ratios between the sizes of the gas disk 
and the stellar disk. We return to this issue in Section 0- 



3.4 The influence of bulge formation 

As already eluded to by Bullock et al. (2000), if the ex- 
cess gas mass in the centres of the disks is transformed into 
a bulge component, this may solve the problem with the 
overly concentrated disks. We now examine whether the in- 
clusion of our particular model for bulge formation, based 
on the idea that self-gravitating disks transform part of their 
material into a bulge, yields exponential stellar disks. 

In Figure H we plot the 2: = surface densities of 
the same models as in Figure y, except that we now have 
included bulge formation (using the parameters listed in 
Section E^). The low angular momentum systems with 
A = 0.028 have formed massive bulges with Mbuige/Mdisk = 
0.7 — 0.8 (i.e., they are more reminiscent of SOs than spiral 
galaxies). The bulge formation has depleted the cold gas in 
the centre to below Ecrit, suppressing subsequent star for- 
mation and resulting in a stellar disk that is remarkably 
exponential. The same applies to the systems with A = 0.06 
(middle panels), for which Afbuigc/Mdisk ~ 0.07, typical of 
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Except for the model with Mvir(O) = 5 X lQ^h~^ Mq and A = 0.129, all stellar disks are more centrally concentrated than an exponential. 
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late-type Sc spirals. Clearly, the formation of a bulge compo- 
nent out of the low angular momentum disk material seems 
fairly successful in producing exponential stellar disks. In 
the case with A = 0.129, however, hardly any bulge compo- 
nent is produced (Mbuige/Afdisk = 0.001). Consequently, the 
more massive models produce, as in the case without bulge 
formation, stellar disks with a central cusp in excess to the 
outer exponential profile. Depending on the radial range in 
which the exponential is fit, the cusp contains on the order 
of 20 to 30 percent of the mass of the stellar disk. 

Recall that we model the infiow of the gas out of which 
the bulge is formed such that no positive gradient is cre- 
ated in the surface density of the cold gas (Section 2.6). 
This is a fairly ad hoc assumption, especially since it is well 
known that the actual HI distribution of disk galaxies often 
reveal a strong positive gradient reflecting a central hole in 
the HI distributionPJ. However, it turns out that the density 
distributions of the stellar disks are fairly insensitive to the 



T It is currently still unclear to what extent these holes are truly 
devoid of gas, or whether they are actually filled with molecular 
or ionized gas. 



details of the inflow-model. If, for instance, we model the gas 
infiow such that all the gas in the central bin is transformed 
into a bulge component before gas at larger radii is used 
(i.e., we resort to a maximally efficient infiow), the bulge 
dominated model galaxies reveal gas disks with large central 
holes, whereas the stellar disks are only marginally different 
from the standard infiow model. Thus, whereas our modeling 
of the bulge formation is admittedly crude, the main point 
is that despite the fact that we have optimized the cooling 
fiow to produce exponential disks in systems with relatively 
large bulges, we are still left with the problem that the high 
angular momentum, low surface brightness systems, with- 
out a significant bulge component, have stellar disks that 
are more centrally concentrated than an exponential. 



3.5 The influence of feedback 

The energy input into the ISM by supernovae can produce 
galactic winds that can drive mass out of the disk. Feedback 
is therefore a natural process that can reshape the density 
distributions of the disks. In order to solve the problem with 
the excess central densities we need the feedback process to 
be relatively more efficient in the centre than in the outer 
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= 0.06 significant 
profiles that are 
significant bulge 



parts of the disks. There are three eflects that determine 
these relative efficiencies within our simple feedback picture: 
At smaller radii more stars form, and thus more SN energy 
is available to drive an outflow. At the same time, the es- 
cape velocities are higher at smaller radii, suppressing the 
efficiency of outflow. Finally, star formation at smaller radii 
starts earlier, when the total escape velocity of the system 
is still low. We now examine how these effects compete and 
whether or not we can produce bulge-less disks with expo- 
nential surface densities. 

In Figure H we plot, for three different feedback efficien- 
cies, the baryonic mass fractions at 2; = in each radial 
bin that have ended up as bulge mass, disk stars, cold gas, 
or have been ejected. In the case without feedback (upper 
panel) , virtually all the gas in the centre of the galaxy ends 
up in the bulge, while at larger radii, most of the gas is con- 
verted into disk stars. The introduction of a modest amount 
of feedback with ea, — 0.02 (middle panel) suppresses both 
the bulge formation and the formation of disk stars. Most 
remarkably, the fraction of the gas mass that is ejected is 
almost constant with radius. Apparently the three effects 
mentioned above cancel each other out such that the frac- 
tional feedback efficiency is virtually independent of radius. 
In the case with £fb = 0.2 (lower panel) the relative feedback 



efficiency is slightly higher in the centre than at larger radii, 
but the effect is not enough to alleviate the problem with 
the central cusps. We thus conclude that our simple picture 
of supernovae induced feedback does not significantly alter 
the density distribution of the resulting disk galaxies (except 
for an absolute offset), and thus can not solve the problem 
of the overly concentrated disks. 



4 TRUNCATION RADII 

As is evident from Figure H, our models reveal distinct trun- 
cation radii in the density distributions of both the stars 
and the gas. The truncation radius of the cold gas reflects 
the maximum specific angular momentum of the baryonic 
mass that has cooled. The truncation radius of the stars, 
on the other hand, refiects the presence of a star formation 
threshold density. 

In order to make a quantitative comparison with the ob- 
served truncation radii in disk galaxies, we construct three 
samples of 100 model galaxies. The samples only differ in 
the value of the feedback efficiency; all other model param- 
eters are kept fixed at their fiducial values listed in Sec- 
tion 2.8, Masses are drawn randomly from a uniform distri- 
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Figure 7. The ratios Rt/Rd (upper panels) and Rui/Rd (lower panels) as functions of log(Vmax)- Results are plotted for three samples 
of 100 model galaxies each, that only differ in the value of the feedback efficiency as indicated above each column. Open (filled) circles 
correspond to models with bulge-to-disk mass ratio less (more) than 0.2. The horizontal dotted lines indicate the range of observed 
values: Rt/Rd = 3.5 it 1.0 and Rui/Rd = 4.1 it 1.8. See text for a detailed discussion. 



bution with 3 x lO^/i"^ Mq < Mvir(O) < 3 x lO^^/i"^ M@, 
the spin parameter is drawn from the probability distribu- 
tion of equation ( pT[ ) with A = 0.06 and ax = 0.6, and the 
MAH parameter a is drawn from a Gaussian distribution 
with a = 0.6 and a a = 0.15. For each model galaxy we de- 
termine the radii 7?,, Rgas, and Rhi- Here we define R, as 
the maximum radius with non-zero surface density of stars, 
Rgas as the maximum radius with non-zero surface density 
of cold gas, and Rni as the radius at which the surface den- 
sity of the cold gas is equal to 1.3 Mopc"^. Taking account 
of the Helium mass abundance, Rm thus corresponds to the 
radius where the HI surface density is equal to I.OMqpc"^. 
In order to express these truncation radii in terms of the 
scale length of the stellar disk we fit an exponential to the 
7-band surface brightness profiles of the models. Since, as 
we have shown above, not all stellar disks are well-fit by an 
exponential, there is some ambiguity involved in these fits. 
After some experimenting we found adequate results upon 
restricting the fits to the radial interval 0.2i?* < R < 0.9R*, 
which is what we adopt throughout. 



4.1 The extent of stellar disks 

The upper panels of Figures M and H plot R,/Rd as func- 
tion of the maximum of the rotation curve and the central 
7-band surface brightness of the best-fit exponential, respec- 
tively. Observations have shown that stellar disks in spiral 
galaxies reveal truncation radii, i?*, in the range from 2.5 
to 4.5 disk scale lengths (van der Kruit & Searle 1981a,b; 
Romanishin, Strom & Strom 1983; Barteldrees & Dettmar 
1994; Pohlen, Dettmar & Liitticke 2000; de Grijs, Kregel & 
Wesson 2001). This range is indicated by horizontal dotted 
lines. It is important to realize that most studies of stellar 
truncation radii have focussed on disk galaxies with small 
bulges. Therefore, in order to make a fair comparison with 
our models, we plot model galaxies with Mbuige/Afdisk > 0.2 
with separate symbols (filled circles). 

In the case without feedback (left panels) the model 
galaxies with A/buigo/Mdisk < 0.2 (open circles) have stellar 
cut-off radii 2 ^ R^/Rd ^ 4.5, in excellent agreement with 
observations. In addition, the models predict that the more 
bulge-dominated galaxies have R^/ Rd ^ 2. It should be rel- 
atively straightforward to test this prediction with deep pho- 
tometry of a sample of Sa and SO galaxies. 

Increasing the feedback efficiency reduces the average 
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Figure 8. Same as Figure M except tliat now Rt/Rd and Rm/Rd are plotted as functions of the / band central surface brigiitness of 
the best fit exponential. 



bulge-to-disk ratio of the models and produces disks with 
lower central surface brightnesses. In the case where Cfb ~ 
0.2 only two of the model galaxies have Mbuige/Mdisk > 0.2. 
Furthermore, the models now predict that low-mass systems 
should have values of R*/Rd ^ 2.5, significantly lower than 
for the more massive systems. Unfortunately most studies 
of stellar truncation radii have focussed on more massive 
disk galaxies, and a more detailed comparison of our mod- 
els with observations has to await more data, especially on 
dwarf and LSB galaxies. Nevertheless, it is reassuring that 
our models with Mbuige/Mdisk ^ 0.2 have stellar trunca- 
tion radii in good agreement with observations, supporting 
the idea that they originate from a star formation threshold 
density set by Toomre's stability criterion. 



4.2 The extent of gas disks 

Little is known about the truncation radii of the gaseous 
component of disk galaxies, mainly due to the fact that the 
gas is generally probed by the HI. Since at column densities 
below ~ 10^^ cm~^ one expects the gas to be ionized by the 
cosmic background flux of ionizing photons (e.g., Sunyaev 
1969; Maloney 1993), a truncation in the HI does not neces- 
sarily reflect the true outer edge of the gas disk. Therefore, 
one typically expresses the size of the gas disk by the radius 



Rui at which the HI surface density is equal to l.OMopc '^ 
(1.2 X 10^° cm-2). 

The lower panels of Figures M and M plot Rni/Rd 
as function of Vinax and /ig, respectively. Systems with 
Mbuigo/Mdisk > 0.2 typically have Rm/Rd ^ 3. For 
the galaxies with lower bulge-to-disk ratios we flnd that 
Rm/Rd = 2.8 ± 0.5 for Knax > 150 kms"^ and Rm/Rd ^ 3 
for Knax < 150 kms~^. For a sample of 73 dwarf galaxies 
with Knax ^ lOOkms^^ Swaters (1999) found an average 
of Rm/Rd ~ 4.1 ± 1.8 (with Rd measured in the J?-band, 
which we assume to be roughly equal to that in the 7-band). 
This is in good agreement with our models with efb = and 
Efb ~ 0.02. However, the model with ejb ~ 0.2 yields and av- 
erage of Rm/Rd = 2.3 for galaxies with Knax < lOOkms"^, 
which is only marginally consistent with the data. 

In Figure M we plot the ratio Rgas/Rm as function of 
Ho for the model without feedback (left panel). Models with 
efb = 0.02 and 0.2 yield very similar results. Except for low 
mass LSB galaxies with /iq i^ 22 mag arcsec"'^ , we find that 
7?gas — -Rhi, indicating that the physical truncation of the 
cold gas occurs close to the radius where Ehi ~ I.OMqpc"^. 
We can compare these model predictions to data by com- 
puting the ratio -RHi.max/^Hi, where i?Hi,max is defined as 
the maximum radius out to which HI is detected. -Rhi, max 
can thus be considered a lower limit on -Rgas. Results for 
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Figure 9. The ratio iJgas/^Hli between the actual cut-off radius of tlie gas disk (iJgas) and the radius where the HI surface density 
equals 1.0 Mq pc~^ (^Hl)) ^^ function of the /-band central surface brightness. The panel on the left shows the results for a sample 



of 100 model galaxies with eft, = 0.0 (all other model parameters are set to their fiducial values, see Section 2.8). Symbols are as in 
Figures pi and pi Note how, except for a few low mass LSB galaxies, all models have Hgas — Rm't the truncation of the model gas disks 
occurs close to the radius where Spji = 1.0 M© pc~^. This is in clear conflict with the data for 26 disk galaxies with readily available HI 
surface brightness profiles, plotted in the right panel. Here -Rgas is the maximum radius out to which HI (or Ha) is detected, and thus 
corresponds to a lower limit on the actual cut-off radius. Data is taken from the following sources: DDO 154 (Carignan & Beaulicu 1989), 
NGC 247 (Carignan & Puche 1990a), NGC 3198 (Begeman 1989), NGC 3109 (Jobin & Carignan 1990), NGC 5585 (Cote, Carignan & 
Sancisi (1991), NGC 7793 (Carignan & Puche 1990b), and Swaters (1999). If no 7-band surface brightness data is available we used the 
colors of de Jong (1996) for conversion. NGC 253 is a special case. This galaxy in the Sculptor group is the only galaxy in which ionized 
gas has been detected at radii beyond i?Hl (Bland-Hawthorn, Freeman & Quinn 1997). With i?Hl = 8.6 kpc (Puche, Carignan & van 
Gorkum 1991) and H« and [Nil] detected out to 11.6 kpc, this implies a lower limit of Rgns/Rm = 1.35. 



a number of disk galaxies with readily available HI surface 
brightness profiles are shown in the right panel of Figure ^ 
As is evident from Figure H, real galaxies have HI disks that 
extend to at least ~ 1.5i?Hi. At these radii the HI surface 
density typically has fallen to ~ 0.1 Mq pc~^. Recall that 
these maximum radii are lower limits to the actual edge 
of the gas disk. Higher sensitivity HI and/or Ha observa- 
tions are required to determine the true extent of the disks. 
For the moment, however, it is sufficient to realize that 
our models, which predict truncation radii close to where 
Ehi ~ 1.0 M0 pc~'^, are in clear conflict with the data. 



5 A NEW PROBLEM FOR THE FORMATION 
OF DISK GALAXIES 

As shown above our models predict that LSB galaxies reveal 
surface brightness profiles that are more centrally concen- 
trated than an exponential. The majority of observed LSB 
disk galaxies, however, have surface brightness profiles that 
are accurately fit by an exponential down to the very cen- 
tre (McGaugh & Bothun 1994; de Blok, van der Hulst, & 
Bothun 1995; Swaters 1999). An exception are the giant, 
massive LSB disk galaxies, which have surface brightness 
profiles that are more concentrated than an exponential, but 
this refiects the presence of a relatively massive bulge com- 



ponent (Sprayberry et al. 1995; Pickering et al. 1997). Our 
failure to reproduce bulge-less LSB disks with exponential 
surface brightness profiles must thus be considered a failure 
of the models. 

In addition, we have shown that whereas our imple- 
mentation of star formation threshold densities yields stel- 
lar truncation radii in good agreement with observations, 
the models predict truncation radii in the gaseous compo- 
nent at too high surface density. Both of these problems 
indicate that our models predict disk galaxies that are too 
concentrated, and in what follows we shall refer to this as 
the disk concentration problem (hereafter DCP). Below we 
examine whether this DCP might be an artifact of the over- 
simplified nature of our models by discussing how various 
modifications of the ingredients or assumptions of our model 
influence the density distribution of the resulting disks. 

5.1 Star formation, galactic winds, and AGNs 

Although we adhere to empirical relations to model the star 
formation in the disk, one could easily envision an alter- 
native model that results in relatively lower star formation 
efficiencies in the central regions. Although this might pro- 
duce stellar disks that more closely resemble exponentials, it 
does not solve the problem at hand. Changing the star for- 
mation efficiencies only changes the relative ratio of gas to 
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Figure 6. The fractions of tlie baryonic mass in eacli radial bin 
that, at z = 0, have ended up as bulge mass (densely cross- 
hatched regions), disk stars (loosely cross-hatched regions), cold 
gas (diagonally-hatched regions), or that have been ejected (white 
areas). All three panels correspond to the fiducial model galaxy 
with Mvir(O) = 5 X 10"h-i Mq, A = 0.06, a = 0.6, and 7 = 1.0, 
but differ in the feedback efficiency: Efj, = 0.0 (upper panel), 0.02 
(middle panel), and 0.2 (lower panel). Note how the baryonic 
mass fraction that is ejected depends only weakly on radius. 



stars, while leaving the central mass densities intact. There- 
fore, such modifications will not alter the circular velocity 
curves, which, at least for the bulge-less LS B g alaxies, seem 
inconsistent with observations (see Section 5.4 below). 

A more plausible solution involves somehow removing 
the excess gas from the centres of the galaxies, or prevent- 
ing it from reaching the centre in the first place. Feedback 
processes seem the most natur al m eans to accomplish this. 
However, as shown in Section 3.5, our simplistic picture of 



SNe induced galactic winds is unable to solve the problem at 
hand. What is needed is a feedback process that is relatively 
more efficient at the centre than at larger radii. It is not in- 
conceivable that a more sophisticated treatment of stellar 



feedback might actually accomplish that. For instance, the 
presence of active galactic nuclei (AGNs), which have been 
ignored in our models, might provide an additional amount 
of energy input exactly at the centre of the galaxies where 
it is most required. However, numerous studies have shown 
that the mass of the AGN's energy source, the massive black 
hole, is proportional to the mass of the bulge in which the 
AGN is embedded (Kormendy & Richstone 1995; Ferrarese 
& Merritt 2000; Gebhardt et al. 2000). Since the DCP is 
most severe for bulge-less systems, where one thus expects 
no (significant) AGN, it seems unlikely that the inclusion of 
AGN related feedback processes can solve the problem with 
the overly concentrated LSB galaxies. Nevertheless, because 
of our extremely limited understanding of the various feed- 
back processes that are likely to play a role in the formation 
and evolution of (disk) galaxies, it is premature to exclude 
feedback as a mechanism that may solve the problem with 
the central concentration of the LSB model galaxies.. How- 
ever, it is very unlikely that feedback processes can actually 
solve the problem with the truncation radii of the gas com- 
ponent, which manifests itself at radii beyond the stellar 
disk, and feedback as a solution to both aspects of the DCP 
thus seems unlikely. 

5.2 Viscosity 

We have not included viscous transport in our models. Vis- 
cosity is efficient in redistributing the disk's angular momen- 
tum distribution, and thus in modifying the density distribu- 
tion of the disk material. Numerous studies in the past have 
in fact argued that, as long as the viscosity time scale is sim- 
ilar to the star formation time scale, viscous transport has 
the natural tendency to produce exponential disks, indepen- 
dent of the initial density distribution of the gas disk (Lin & 
Pringle 1987; Yoshii & Sommer-Larsen 1989; Clarke 1989; 
Olivier, Blumenthal & Primack 1991; Ferguson & Clarke 
2001). However, there are two important caveats here. First 
of all, in most cases the resulting stellar disk has a density 
distribution that is only exponential in the outer regions, 
while the central regions reveal a central cusp, reminiscent 
of our disk profiles. Secondly, virtually all studies of viscous 
disks start with an initial density distribution that is less 
concentrated than an exponential. However, we have shown 
that the standard model for disk formation predicts density 
distributions that are already too concentrated, even before 
the onset of viscosity. Finally, it is worth stressing that the 
main mechanism for viscosity is poorly constrained, and in 
particular, that there is no detailed theory that links the 
time scale of viscosity to that of star formation, as required. 
Most importantly, viscous transport is oriented inward 
in the central regions and outward in the outer regions. 
Thus, whereas viscosity seems the natural solution for the 
problem with the truncation radii, it will at the same time 
only aggravate the problem with the overly concentrated 
LSB galaxies. We thus do not consider viscosity a viable 
solution to the DCP. 



5.3 The mass accretion histories of dark matter 
haloes 

The DCP outlined above only applies to LSB galaxies. 
Throughout, it has been assumed that LSB galaxies only 
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Figure 10. The influence of the MAH on the stellar surface den- 
sity profiles of the disks. Plotted are log(Estar(R)) for four model 
galaxies with Mvir(O) = 5 X 10^^ h'^ Mq and A = 0.129. The 
models only differ in the values of Zf and a that parameterize 
the MAH. Note how despite a wide range in MAHs, the resulting 
disks all are too centrally concentrated to be consistent with ob- 
servations. This illustrates that the DCP is not solved by resorting 
to specifically tuned mass accretion histories. 



distinguish themselves from HSB galaxies in a larger spin 
parameter. Several studies have shown that this assumption 
is consistent with a wide variety of observed properties of 
LSB galaxies (e.g., Dalcanton et al. 1997; Mo et al. 1998; 
Jimenez et al. 1998). One important shortcoming of these 
models, however, is that they offer no explanation for the 
fact that LSB galaxies are less strongly clustered than their 
HSB counterparts (Mo, McGaugh & Bothun 1994). This has 
prompted several authors to suggest that differences in disk 
surface density are driven by differences in the amplitude 
of the original density fluctuations from which the galaxies 
arise (e.g., McGaugh & de Blok 1998). In this case, HSB 
and LSB galaxies are expected to have correlation functions 
with different amplitudes. 

If indeed peak height is the main parameter that dis- 
criminates LSB from HSB galaxies one expects the two to 
have different MAHs (the amplitude of density perturba- 
tions is directly related to the halo formation epoch). It is 
therefore worthwhile to examine whether the disk concentra- 
tion problem for LSB galaxies may be solved by resortin g to 



different MAHs. Furthermore, as emphasized in Section 2.2 



we have limited ourselves to a fairly small subset of possible 
MAHs by only considering models with Zf = 10.0 and a in 
the range [0.3,0.9]. Clearly, a more extended exploration of 
parameter space is needed before we can actually claim a 
failure of the models to produce realistic LSB galaxies. 

In Figure Fui we plot the stellar surface densities of four 
models with Mvir(O) = 5 x 10"/i~^ Mq and A = 0.128 (all 
othe r pa rameters are set to the fiducial values listed in Sec- 
tion 2.8). The models only differ in their values of Zf and a 
as indicated. As is evident, the MAH mainly influences the 



the central density of the disk. All four models, which have 
wildly different MAHs, yield disks that are inconsistent with 
the observed surface brightness profiles of LSB disk galax- 
ies. In particular, the case with Zf = 1.0 (dotted line), which 
is the most appropriate for a scenario in which LSB galax- 
ies form from low amplitude density fluctuations, yields a 
surface brightness profile that is only approximately expo- 
nential over a very small radial range. Clearly, the DCP can 
not be solved by resorting to specifically tuned MAHs. Note, 
however, that in principle haloes that form later will also be 
of lower density, something not taking into account in these 
models (see discussion in Section 2.2). We address this in- 
fluence of halo density on the surface brightness profiles of 
the disks in the next section. 



5.4 The nature of dark matter 

In the standard picture for disk formation, outlined above, 
the density distribution of the disk is related to both the den- 
sity and angular momentum distribution of the dark matter 
halo. If the dark matter is cold and coUisionless it virial- 
izes to produce strongly concentrated haloes with steep cen- 
tral cusps. The disks that form in these haloes will there- 
fore also be strongly concentrated. A possible solution to 
the DCP, therefore, might be that dark matter haloes are 
less concentrated, e.g., they have a constant density core. 
It is noteworthy in this context that numerous studies in 
recent years have claimed that the observed rotation curves 
of dwarf and LSB galaxies imply dark matter haloes with 
constant density cores (Flores & Primack 1994; Moore 1994; 
McGaugh & de Blok 1998, but see van den Bosch et al. 2000 
and van den Bosch & Swaters 2001). Furthermore, if LSB 
galaxies for m from low amplitude density fluctuation (see 
Section 5^), one expects LSB galaxies to be embedded in 
low-density dark matter haloes Therefore, it is worthwhile 
to explore to what extent a modiflcation of the dark matter 
density distribution can solve the DCP. 

The two parameters in our models that set the dark 
matter density distributions are the central cusp slope 7 
and the normalization of the halo concentrations Co. In the 
upper left panel of Figure hi] we plot the stellar surface 
densities of four models with Afvir(O) = 5 x 10^^ h~^ Mq, 
X = 0.128 and a = 0.6. Our standard ACDM model with 
7 — 1.0 and Co = 4.0 (solid lines) yields a nearly bulge-free 
disk (Afbuigc/Mdisk — 0.001) with a central cusp (cf. upper 
right panel of Figure ™. The disk is so strongly concen- 
trated that the rotation curve (upper right panel) reveals 
a steep rise followed by a modest decline. This is another 
manifestation of the DCP, as the observed rotation curves 
of LSBs reveal a slow rise in the central parts and virtually 
never a declining part. The dot-dashed line shows the sur- 
face brightness profile for the same model galaxy but with 
Co = 1.0. This results in a halo concentration parameter 
c that is a factor four times smaller (c = 3.7 compared to 
c = 14.8). Eventhough the mass of the resulting galaxy is 
significantly less concentrated, as is evident from the rota- 
tion curves, the surface brightness of the resulting stellar 
disk is only marginally less concentrated. The main effect 
of reducing the halo concentration is to produce a more ex- 
tended disk, and to increase the bulge-to-disk ratio (lower 
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Figure 11. Properties of model galaxies witli Mvir(O) = 5 X lQ^^h~^ Mq, A = 0.129 and a = 0.6 for four different dark matter density 
distributions. The solid lines correspond to the standard ACDM haloes, with 7 = 1.0 (i.e., haloes have NFW density profiles) and 
Co = 4.0. The dot-dashed lines correspond to a model with a halo concentration parameter that is four times smaller (Co = 1.0). Dashed 
lines correspond to a model with a constant density core (7 = 0.0), but with the same concentration parameter as our fiducial ACDM 
model. The dotted lines, finally, corresponds to a constant density cores with a core radius four times larger than for the standard model. 
The upper panels plot the surface densities of the stellar disks and the circular velocity curves (out to the truncation radius of the cold 
gas) as function of radius. The lower left panel plots the bulge-to-disk ratios (in mass) as function of redshift. Note how less concentrated 
dark matter haloes yield more extended disks with close-to-exponential density distributions and rotation curves that rise more slowly, 
in better agreement with observations. However, because the disks are more self-gravitating they also produce larger bulges, and the 
problem of producing bulgc-lcss LSB disk galaxies with exponential stellar disks remains. Finally, the lower right panel plots the ratio 
Rgms/Rm as function of /Xq for 100 model galaxies with 7 = 0.0 and Co = 1.0. Symbols are as in Figure bl Note how even for dark 
matter haloes with large constant density cores the truncation radii of the gas disks in HSB systems occurs close to i?Hli in conflict with 
observations. 



left panel) . Clearly, merely changing the halo concentrations 
is not able to solve the DCP. 

The dashed lines in Figure hll show the results for the 
same model galaxy but with Co = 4.0 and 7 — 0.0, i.e., 
rather than changing the halo concentration parameter, we 
now consider dark matter haloes with a constant density 
core rather than a r~^ cusp. The results of lowering 7 are 
fairly similar to those of lowering Co: the resulting stellar 
disk is somewhat less concentrated and more extended than 
for the case with the NFW dark matter halo, but the stellar 



density distribution is still significantly more concentrated 
than a pure exponential. Furthermore, the bulge-to-disk ra- 
tio is an order of magnitude larger than in the standard case 
with 7 = 1.0 and Co = 4.0. 

The dotted lines, finally, correspond to a model with 
7 — 0.0 and Co = 1.0. In this case the stellar disk is close 
to exponential over a relatively large radial range, and the 
central cusp has largely disappeared. Furthermore, the re- 
sulting rotation curve is still rising at the outer edge of the 
gas disk, in better agreement with observations. However, 
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as can be seen from the lower left panel, the bulge-to-disk 
ratio is, with Afbuigo/A/disk ~ 0.1, two orders of magnitude 
larger than for the ACDM model with 7 = 1.0 and Co = 4.0. 
This emphasizes the robustness of the DCP. When trying to 
lower the central densities of the disks by resorting to dark 
matter haloes of lower (central) densities, the disks become 
more self-gravitating, which results in larger disk-to-bulge 
ratios. The problem of producing bulge-less exponential stel- 
lar disks therefore remains 

Finally, the lower right panel plots Rga.s/RHi as func- 
tion of Ha for a sample of 100 model galaxies with 7 = 0.0 
and Co = 1.0. Compared to the ACDM model with 7 = 1.0 
and Co = 4.0 (left panel of Figure H), both the fraction 
of low surface brightness galaxies and that of systems with 
Af bulge /Ajfdisk > 0.2 have increased. However, systems with 
fj,Q ^ 22 mag arcsec~^ still have Rgas/Rni — 1, in disagree- 
ment with observations. Modifying the structure of the dark 
matter haloes can thus not solve the problem with the trun- 
cation radii outlined in Section 14 



5.5 The angular momentum 

Under the assumption of angular momentum conservation, 
the density distribution of the disk that forms is a di- 
rect reflection of the angular momentum distribution of the 
baryons in the protogalaxy. The DCP might therefore be 
a consequence of our specific treatment of the evolution of 
the angular momentum content. In particular, we have made 
four crucial assumptions: (i) the specific angular momentum 
of the baryons is conserved, (ii) the spin parameter of the 
dark matter haloes is constant with time, (iii) the angular 
momentum vectors of all mass shells are aligned, and (iv) the 
baryons and dark matter have the same density and angular 
momentum distribution 

Relaxing assumptions (i) and/or (iii) in general will 
only worsen the DCP. As discussed in Section ^ conserva- 
tion of specific angular momentum is required to produce 
disks with the correct size distribution. Furthermore, gas 
with misaligned angular momentum vectors has to align it- 
self in order to form a disk, which it can only do by transfer- 
ring angular momentum to the halo, thus resulting in even 
more compact disks. 

Since disks form from the inside out, one obvious way of 
preventing too much mass from cooling to small radii is to 
relax (ii) and assume that the spin parameters of protogalax- 
ies were systematically higher at higher redshift. However, 
although for each individual halo A is likely to vary with 
redshift, such a systematic trend with z seems to be ruled 
out by the fact that A^-body simulations have shown that the 
spin parameter distribution is virtually independent of mass, 
environment and/or redshift (Lemson & Kauffmann 1999). 
Analytical studies by Heavens & Peacock (1988) and Cate- 
lan & Theuns (1996) have suggested a small anti-correlation 
between the average spin parameter and the peak height of 
density fluctuations. However, the amplitude of this correla- 
tion is much smaller than the spread in A at any given value 
of the peak height, and this effect will thus not significantly 
affect our results. 

We can further check the validity of assumption (ii) by 
comparing the resulting angular momentum distributions of 
the dark matter haloes to A'^-body simulations of structure 
formation. In a recent paper, Bullock et al. (2000, hereafter 
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Figure 12. The mass distribution of specific angular momentum 
for models with Mvir(O) = 5 X 10"/i-l Mq and A = 0.06. Results 
are shown for three different MAHs: a = 0.3 (crosses), a = 0.6 
(open circles), and a = 0.9 (open triangles). The solid lines are 
the best fit profiles of the form (pl|). Shown for comparison is 
the distribution for a uniform sphere in solid body rotation (thin 
dotted curve). The best fit values for fi are listed in the figure, 
and vary from fi = 1.6 (a = 0.9) to fi = 1.9 (a = 0.3), well 
in the range of /j. values found by BOO. This implies that our 
assumptions regarding the evolution of the angular momentum 
result in realistic angular momentum profiles. 



BOO) have shown that dark matter haloes in a high resolu- 
tion ACDM A'^-body simulation have specific angular mo- 
mentum profiles that are well fit by 



^U) 



M{< j) 

Mvir 



(j/jmax) 



(31) 



Here M{< j) is the halo mass with specific angular mo- 
mentum less than j, Mvir is the halo's virial mass, /x is a 
free parameter, and jmax is the maximum specific angular 
momentum in the halo. BOO have shown that the angular 
momentum content of a dark matter halo is well described 
by the pair (A,/i), and that for 90 percent of the haloes 
1.06 < jj, < 2.0 with a mean of (h) = 1.25. We have com- 
puted m{j) for several of our models. The results for three 
models with Mvir(O) = 5 x 10^^h~^ Mq are shown in Fig- 
ure I12I together with the best fit profiles of the form ( pl[ ) . 
As can be seen, the mass distribution of specific angular 
momentum in our models is well described by equation ( |3l[ ) 
with best fit values for h in the range 1.6 ^ ^^ ^ 1.9, in 
good agreement with the results of BOO. We have also tested 
whether h is correlated with mass, and find a very weak de- 
crease of fi with halo massn, again in good agreement with 
BOO. 

Although our best-fit values of /i are well inside the 
range of /i values found by BOO, the spread we obtain is 

T for a = 0.6 we find that fj, decreases from 1.76 for Mvir(O) = 
5 X lO^h-^ M0 to /i = 1.68 for Mvir(O) = 5 x lO^^/j-i j^q 



© 2000 RAS, MNRAS 000, 



18 Frank C. van den Bosch 



0.8 - 



S 0.6 - 



0.4 



0.2 







1 1 1 


/^'-'-^^ ; 


— // 


_ 


/ / 


- 


/ / 


- 


- / / 

- / / 

; / / 

- / / 
-/ / 


- 


- 


- 


'J 


At=106 ■ 


H 


/J - 1 7n 




1 1 1 


1 , , , 1 , , , 1 , " 







4 

j/Jlol 



6 



Figure 13. Mass fractions Af (j)/Mvir(0) as function oi j /jtot for 
two different values of fj,. Decreasing /i increases jmax, and thus 
the amount of high angular momentum material, and simultane- 
ously increases the amount of low angular momentum material. 
Lowering ^ thus has the same effect as viscosity. 



much smaller. However, this is most likely a reflection of the 
fact that our model haloes have smooth, spherical density 
distributions without any substructure. This also explains 
why BOO find a weak correlation between A and /i, whereas 
in our models the m[j) profiles are independent of the value 
of the halo spin parameter. Nevertheless, it is worthwhile 
to examine whether the DCP can be solved by resorting to 
haloes with other values of /i. In order to address this issue 
it is useful to rewrite equation (pll) in the form 



m{j) = fj, 



ij/jt. 



ij/jtot) + (m - l)(imax/itot) 



Here 



Jtot — 



(Vvir 



= Jn 



1-M 



1 - (^ - l)ln 



M-1 



(32) 



(33) 



is the total specific angular momentum of the halo. Sub- 
stituting equation ( |33[ ) in (p2|) one can compute m{j) as 
function of j/jtot for a given value of ^. Results for /x = 1.7 
(the average value found for our haloes) and fi — 1.06 (the 
5 percent lower bound of the distribution found by BOO) are 
shown in Figure nsl As is immediately evident, lower val- 
ues of /i result in higher values of jmax, and thus in more 
extended disks. However, at the same time, decreasing /i in- 
creases m{j) at low values oi j/jtot, which implies that the 
disk will become more centrally concentrated, thus aggra- 
vating the DCP. In fact, lowering /i mimics the effects of 
viscosity (cf. Zhang & Wyse 2000). 

It thus becomes clear that in order to solve the DCP 
it is required that the baryons have an angular momentum 
distribution that is completely decoupled from that of the 
dark matter (i.e., assumption (iv) is in error). The baryons 
need an angular momentum distribution with both less low 



angular momentum material and more high angular momen- 
tum material (i.e., a tail that extends to higher jmax). From 
the above it is clear that this can not be realized by merely 
changing /i. Rather, the baryons need an angular momen- 
tum distribution that is completely different from that of 
equation pl| (see also van den Bosch, Burkert & Swaters 
2001). 



6 DISCUSSION & CONCLUSIONS 

In our current paradigm of galaxy formation, the angular 
momentum of protogalaxies originates from tidal torques of 
nearby density perturbations. It has become clear that this 
mechanism provides just enough angular momentum to pro- 
duce disk galaxies of the right size (Fall & Efstathiou 1980; 
Dalcanton et al. 1997; Mo et al. 19998), implying that the 
baryons can not lose a significant fraction of their specific 
angular momentum. In the standard picture of disk forma- 
tion it is therefore assumed that the baryons conserve their 
specific angular momentum when cooling to form the disk. 
Under these conditions one can directly compute the density 
distribution of the resulting disk once the actual distribution 
of specific angular momentum of the protogalaxy is known. 

Recently, Bullock et al. (2000) have determined the spe- 
cific angular momentum distribution of a large set of dark 
matter haloes in a ACDM universe, and have shown that 
if baryons have the same angular momentum distribution 
as dark matter, the resulting disk galaxies are more cen- 
trally concentrated than an exponential. This then raises 
the question whether or not baryonic processes related to 
star formation, bulge formation, and feedback can produce 
stellar disks with exponential surface brightness distribu- 
tions out of such highly concentrated gas disks. In this pa- 
per we have presented new models for the formation of disk 
galaxies to address this question. We use a simple parame- 
terized description of the mass accretion histories of the dark 
matter haloes, and, following Firmani & Avila-Reese (2000) 
and Avila-Reese & Firmani (2000), the assumption is made 
that the halo spin parameter is constant with time. We have 
shown that this implies specific angular momentum distri- 
butions that are in good agreement with those obtained by 
Bullock et al. (2000) from high resolution A'^-body simula- 
tions. Including cooling and adiabatic contraction, we have 
confirmed the results of Dalcanton et al. (1997), Firmani & 
Avila-Reese (2000), and Bullock et al. (2000) that the re- 
sulting gas disks are more centrally concentrated than an 
exponential. Next we included simple prescriptions to de- 
scribe star formation, bulge formation, feedback, and chem- 
ical evolution, and investigated whether these processes can 
transform the highly concentrated gas disks into stellar disks 
with close to exponential surface brightness distributions as 
observed. 

At first sight our models are remarkably successful in 
producing disk galaxies with density distributions as ob- 
served. The essential ingredient is bulge formation, which 
prevents the formation of systems with strongly declining 
rotation curves. Furthermore, the process of bulge forma- 
tion results in stellar disks with close to exponential surface 
brightness profiles. This result is fairly insensitive to the de- 
tails of the bulge formation process, as long as it is coupled 
to a stability criterion for the disk (cf. van den Bosch 1998). 
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The introduction of a star formation threshold density, as 
motivated by the results of Kennicutt(1989), yields gas mass 
fractions and stellar truncation radii in excellent agreement 
with observations (cf. van den Bosch 2000; van den Bosch 
& Dalcanton 2000). 

Despite these clear successes of the model, a closer in- 
spection of the model galaxies reveals two important short- 
comings. The first problem mainly concerns LSB disk galax- 
ies. These systems form out of protogalaxies with high an- 
gular momentum and do not produce a significant bulge 
component. This fact that LSB galaxies have typically lower 
bulge-to-disk ratios than their HSB counterparts is in good 
agreement with observations. However, the resulting stellar 
disks of these LSB systems are too centrally concentrated. 
Although their surface density profiles are close to expo- 
nential at the outside, they reveal a strong central cusp in 
clear disagreement with observations. The second problem 
concerns the extent of the gas disks. The models predict 
truncation radii of the gas close to the radii Rui where the 
HI column density has fallen to ~ 10^" cm^^. In real disk 
galaxies, however, HI is observed out to radii well in excess 
of -Rhi, again in clear contradiction with the models. 

The problem with the truncation radii seems straight- 
forward to solve by including viscosity, which will transport 
disk mass in the outer parts to larger radii. However, at 
the same time viscous transport is oriented inwards at small 
radii. Thus, whereas viscosity seems the obvious solution for 
the problem with the truncation radii, it will only aggravate 
the problem with the central concentration of (LSB) disks, 
and additional processes are required to solve the problems 
at hand. 

Including a simple model for a galactic wind induced 
by SNe can expel large amounts of baryonic matter from 
the disk, but it does so with a relative efficiency that is vir- 
tually independent of radius. This implies that the actual 
density distribution is, except for an offset, left intact. Ad- 
ditional energy input from AGNs, not taken into account in 
our models, is unlikely to solve the problem with the overly 
concentrated LSB galaxies, since the mass of massive black 
holes is strongly correlated with that of the bulges, which 
are virtually absent in LSB galaxies. 

Some studies in the past have suggested that LSB galax- 
ies differ from their HSB counterparts in that they form in 
density peaks of lower amplitude, rather than in peaks with 
more angular momentum (as assumed here). This implies 
that LSB galaxies form later and inside haloes of lower den- 
sities. By exploring a wide range of MAHs and halo concen- 
trations we have shown that even for this picture the LSB 
disks are too centrally concentrated. 

We therefore conclude that understanding the forma- 
tion of disk galaxies in CDM cosmologies faces two impor- 
tant challenges. We first of all need a mechanism that can 
prevent the angular momentum catastrophe which results in 
disks being an order of magnitude too small. Furthermore, 
even if the mass accretion is smooth and the angular mo- 
mentum conserved, the disks that form, although of the right 
size, are too centrally concentrated. The robustness of this 
problem seems to suggest that the baryons need an angular 
momentum distribution that is clearly distinct from that of 
the dark matter. Similar results were recently obtained by 
Navarro & Steinmetz (2000) and van den Bosch, Burkert 
& Swaters (2001). Navarro & Steinmetz used simple scaling 



relations to show that, in a ACDM Universe, disk galaxies 
only need to accrete a small fraction of the total baryonic 
mass to match the zero-point of the TuUy-Fisher relation, 
but must draw a comparably much larger fraction of the 
available angular momentum. Van den Bosch, Burkert & 
Swaters computed the angular momentum distributions of 
a sample of low-mass disk galaxies from the observed ro- 
tation curves and disk density distributions. A comparison 
with the angular momentum distributions of cold dark mat- 
ter haloes found by Bullock et al. (2000) , clearly reveals that 
disks lack predominantly low angular momentum material 
compared to their dark matter haloes. 

An interesting alternative to considering a decoupling 
between the dark and baryonic mass components, is to 
change the nature of the dark matter. Recently, numer- 
ous studies have focussed on scenarios in which the dark 
matter is warm (WDM) or self-interacting (SIDM) . In both 
cases one expects dark matter haloes to have constant den- 
sity cores with less substructure than in CDM models (e.g., 
Spergel & Steinhardt 2000; Bode, Ostriker & Turok 2001). 
This not only alleviates the problems with disk rotation 
curves and the angular momentum catastrophe, but it will 
also result in disks that are less centrally concentrated. In 
order to investigate whether a modification of the nature 
of dark matter might solve the disk concentration prob- 
lem, we have investigated models in which the haloes have 
(large) constant density cores. As expected, the resulting 
disks are more extended as in the CDM case, and with sur- 
face brightness profiles that start to approach an exponen- 
tial form. However, at the same time the disks become more 
self-gravitating, resulting in the formation of relatively mas- 
sive bulges. Therefore, the problem of producing bulge-less 
LSB galaxies with exponential surface brightness profiles re- 
mains. Before concluding that the disk concentration prob- 
lem persists even in WDM and SIDM scenarios, we caution 
that throughout we have made the assumption that the halo 
spin parameter does not evolve with time. Although this 
seems a valid assumption to make in the case of CDM, very 
little is known about the (evolution) of the angular momen- 
tum distribution of haloes in alternative dark matter scenar- 
ios. An investigation of the distribution of angular momen- 
tum in haloes of warm and/or self-interacting dark matter 
will proof extremely useful in this context. 
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